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ABSTRACT 

We report the detection of the CO J = 1 — emission Hne in [H89] 1821+643, one of the most 
optically luminous QSOs in the local Universe and a template ULIRG-to-QSO transition object, 
located in a rich cool-core cluster at 2; = 0.297. The CO emission is likely to be extended, highly 
asymmetric with respect to the center of the host elliptical where the QSO resides, and corresponds 
to a molecular gas mass of ^ 8.0 x 10^ M©. The dynamical mass enclosed by the CO line emission 
region could amount up to ~ 1.7 x 10^^ Mq (80% of the total mass of the elliptical host). The bulk 
of the CO emission is located at ~ 9 kpc south-east from the nuclei position, and close to a faint 
optical structure, suggesting that the CO emission could either represent a gas-rich companion galaxy 
merging with the elliptical host or a tail-like structure reminiscent of a previous interaction. We argue 
that the first scenario is more likely given the large masses implied for the CO source, which would 
imply a highly asymmetric elliptical host. The close alignment between the CO emission major axis 
and the radio-plume suggests a possible role of the latter excitation of the ambient gas reservoir. The 
stacking technique was used to search for CO emission and 3mm continuum emission from galaxies 
in the surrounding cluster, however no detection was found either toward individual galaxies or the 
stacked ensemble of galaxies, with a Scr limit of < 1.1 x 10^ Mq for the molecular gas. 
Subject headings: galaxies clusters: general — galaxies:formation — galaxies: evolution — galaxies 

individual (J1821+643) — galaxies: ISM — galaxies:starburst — galaxies: QSOs: 

general 



1. INTRODUCTION 

In the general picture for QSO and stellar spheroid 
formation, two gas-rich disk galaxies merge and the 
large amount of gas, compressed and shocked by the in- 
teraction, trigger intense starburst activity in the sys- 
tem, which is seen as a powerful IR luminous galaxy 
(so called Ultr a Luminous IR Galaxies ~ ULIRGs; 
iSanders fc Mira bel 1996). As the system evolves, the gas 
flows into its central region, likely feeding an active galac- 
tic nuclei (AGN). When the black- hole has grown enough 
in size and luminosity, feedback and gas consumption 
from star-formation and AGN activity eventually dissi- 
pates the dust an d an optically luminous QSO emerges 
(|Di Matteo et al.l 120051 ). Once AGN activity subsides, 
the galaxy evolves into a massive spheroid, hosting a 
quiescent super-massive black-hole (SMBH) at its cen- 
ter. In this scenario, optically and IR luminous QSOs 
with warm IR colors (S'25^m/<5'60Mm > 0.2) represent the 
rare U LIRG-to-QSO transition objects (jSanders et al.l 
Il988allbh . 

Observations of CO line emission found large amounts 
of molecular gas associated with these ULIRG-to- 
QSO systems , in support o f an evo lutionar y sce- 
nario (e.g. Sanders et aLl 119891: lYoung fc Scoville^ 119911: 
ISanders fc Mirabell 11996^ . However, due to the neces- 
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sarily low resolution of single-dish telescopes, poor sensi- 
tivity of the few interferometric CO maps available and 
the lack of high-resolution optical images (e.g. with the 
HST) for most of these objects, few studies were able 
to determine the location of t he molecular g a s rese rvoir 
respect to the powerful QSO. IScoville et al.l (|2003D im- 
aged the CO line emission from a number of the most 
optically luminous PC QSOs, which are thought to be 
the end-product of the starburst-to-QSO transition. In- 
terestingly, large amounts of gas, similar to that found 
in ULIRGs, where found in these objects. The conclu- 
sion drawn from these observations was that, assuming 
that the CO emission and starburst activity are well cen- 
tered around the nuclei, the QSO ho st galaxies should be 
late- type spirals (jScoville et al.ll2003l ). However, this con- 
tradicts the fact that most of these objects are typically 
found to have giant elliptical host galaxies in the optical. 

Recently, sensitive high-resolution interferometric CO 
observations of a template ULIRG-to-QSO transition 
system, HEG450-2958 at z 0.3, along with Hubble 
Space Telescope (HST) imaging found a merger between 
a gas-rich spiral and a gas-poor elliptical, instead of an 
optically luminous Q SO hosted by a late-type gas-rich 
spiral in this system (jPapadopoulos et"aril2008a|) . Here, 
the optically luminous QSO lies outside the molecular gas 
reservoir, which is rather ide ntified with a massive but 
obscured compa nion galaxy (Pa padopoulos et"aLll2008al : 
lElbaz et "aLl 120093 . High-resolution IR observations of 
this object support this assertion, and suggest also a 
more complex scenario where the QSO radio jet alters the 
state of the ISM and possibly induc es a starburst in the 
gas-rich member of the interaction (|Klamer "eral]l200l 
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TABLE 1 

Observed properties for the H1821+643 system 



Source 


"CO 


■500 


z " 






Mgas " 


SgsGHz 




(J2000) 


{J2000) 




(mjy beam ^) 


(Jy km s"-*-) 


(XIO^ Mq) 


(mJy) 


H1821+643 






0.2970 


< 2.1 


< 1.0 




10.3 ±0.5 


Companion t 


. . . 18 21 57.42 


+64 20 34.75 


0.2974 


2.1 ±0.5 


2.3 ±0.4 


(8.0 ± 1.7) 


< 1.5 



t CO emitting source. 

" redshift for H1821+643 nuclei taken from HSchneider et al.|[T99l ) 
* Peak emission measured in a CO image averaged over the line FWHM. 
Velocity and area integrated CO emission. 

Gas mass computed from the measured CO integrated flux and using a Xco prescription as explained in the text. 
Continuum flux density at 93 GHz. 



iFeain et all 120071: lElbaz et all I2OO90 . Gas-poor /gas-rich 
galaxy interactions "igniting" powerful QSOs and star- 
bursts are expected to be common under the hierarchica l 
framework of structure formation ('S pringel et al.|[2005D . 
Moreover, an important consequence of unrelaxed dy- 
namical merger configurations between a gas-rich system 
and a large QSO-hosting elliptical would be the serious 
under-estimation of the dynamical mass of the latter, 
an effect already noted for more orderly ( but still un- 
relax ed) gas-rich spirals at high redshifts (jPaddi et al.l 
I20T0I) . This can have direct impact on CO-based es- 
timates of Q SO host galaxy masses at high redshifts 
(|Wang et al.ll2010D if they involve such mergers. Thus, 
such reported ULIRG-to-QSO transition objects must be 
revisited with high resolution CO and optical imaging, in 
order to reveal the detailed dynamic configurations of the 
molecular gas reservoirs with respect to the AGN, and in- 
vestigate the possible role of the latter in the formation 
of galaxies. For local powerful QSOs, typically hosted 
by well-characteriz ed large ellipticals with known masses 
(jFlovd et al.ll2004l ). high resolution CO observations can 
serve as excellent test-beds of the accuracy of dynamical 
mass estimates using CO lines. 

One of the most prominent local templates of candi- 
date ULIRG-to-QSO transition objects is IR/optically 
luminous QSO [HB89] 1821+643 system (hereafter 
H1821+643) at z = 0.297. In this paper, we report CO 
emission line observations of this system using the Com- 
bined Array for Research in Millimeter-wave Astronomy 
(CARMA). These observations directly trace the bulk of 
molecular gas, which is the fuel for star-formation. We 
use a concordance ACDM cosmology throughout, with 
ffn = 71 km s"^ M pc-\ VL^ = 0.27 and VLt^ = 0.73 
HSpergel et al.l[200l . 

2. OBSERVATIONS 

2.1. H1821+643 

With an optical absolute magnitude My = —27.1 at 
a redshift of 0.297, H1821+643 is one of the most lumi- 
nous QSOs in the loc al universe (IHutchines fc NefalTMl: 
[Schneider et al.lll992[ ). H1821+643 has been subject to a 
number of studies aimed at determining the morphology 
of its host galaxy. This bright QSO is hosted by a giant 
elliptical galaxy, with a half-light radius of ^ 14 kpc in 
extent, as revealed by modeling of th e point-spread func- 
tion ( PSF) subtra cted HST images (iMcLeod fc McLeodI 
120011: iFlovd et al.l I2004D . The luminous nuclear com- 
ponent outshines the underlying host with a luminos - 
ity ratio of - 11 in the HST I-band (,Flovd et al.|[200l . 



111821+643 lies i n a massive galax y cluster with optical 
richness class 2 (|Lacv et al.l Il992i) . X-ray observations 
show this corresponds to a strong cool-core cluster with 
a cooling time of 1 Gyr, a decrease of temperature in 
the cluster center from 9.0 ± 1.5 to 1.3 ± 0.2 keV, and 
a mass a ccretion rate fro m the central engine of 40 Mq 
yr"i (Russell et al.l[2010l) . 

Radio continuum imaging of this 111821+643 indicates 
this is a radio-quiet QSO, with a radio-lobe extending 
out t o 250 kpc from the bright c ore in south-west direc- 
tion (Papadopoulos et al.' 1995; Blundell & Lacy' 1 19951 : 
Blundcll et al. 1996; Blundell & Rawhngs 2001). The 
radio-lobe appears to bend by ^ 80° in projection on sub- 
arcsec scales, suggesting that the radio axis is precess- 
ing, which implies a SMBH binary system at its nuclei 
(|Blunde ll fc Rawlings' '2001") that has possibly already 
merged (Robinson et al. 2010). Interestingly, this object 
has radio luminosities at the transition to separate ra- 
dio quiet and radio loud QSOs and also at the transition 
betw een FRI and FRII structures (jBlundell fc Rawlingsl 
I2OOI . 

2.2. CO observations 

Observations of the ^^CO J = 1 — emission line 
(rest frequency: 115.271 GHz) were performed using 
CARMA in C-configuration with two overlapping tar- 
get fields. The first field was observed during October 
to November 2008 centered at aj2ooo = 18''21'"57.21^ 
'5j20oo — 64°20'36.23" (6 tracks; centered at the posi- 
tion of H1821+643). Observations of the second field 
were obtained during March 2010 in 2 tracks centered 
at a = 18''22"02.8^ S = 64°20'05.3", - 47.7" south- 
east of the first pointing. We offset the second point- 
ing to allow us to simultaneously observe H1821+643 
plus another source of interest in the field. This object 
will be studied in a separate publication. The observa- 
tions were taken under good millimeter weather condi- 
tions with r230GHz < 0.5. 

We used the 3 mm receivers tuned to 88.8751 GHz 
(lower sideband). Three bands (or spectral windows) of 
15 channels each and 31.25 MHz per channel were used. 
These bands were put together to roughly cover ^ 1.3 
GHz in frequency. For the 2010 observing runs, the 
three bands were overlapped by 2 channels in order to 
improve the bandpass calibration, and the edge channels 
were flagged. In the 2008 runs, however, the bands were 
not overlapped and thus no flagging of the edge chan- 
nels could be performed as it would translate in spectral 
gaps. This resulted in a loss of sensitivity in channels 
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Fig. 1. — Left: The HST I-band image is shown in the background with green and red contours overlaid showing the 3mm continuum 
and the CO 1-0 map integrated over 540 km s~^, respectively, towards H1821+643. Green contours are shown in steps of 8(7, with u = 0.2 
mjy beam~^, starting at ±8cr. Red contours levels are shown in steps of la = 0.5mjy beam^^, starti ng at ±2a. The and "x" white 
symbols represent the position of the bright radio nuclei and south-west radio component identified bv'B lundell et ahl 119961) . The dashed 
orange line represents roughly the direction of the jet (from Blundell et al. 1996; Blundell & Rawlings 2001]). The red ellipse at the bottom 
left shows the beam shape and size. Right: CO 1-0 spectrum of H1821-t-643 averaged over a region that encloses the CO emission shown 
in the left panel. The solid line represents a Gaussian fit to the data.. 



with velocities of ~ ±780 km s^^. After combining all 
datasets, the effective bandwidth was ~ 1.281GHz. The 
large bandwidth allowed us to search for emission not 
only from H1821-f 643 but also from galaxies in its dense 
environment. The upper sideband was tuned to ^ 93.1 
GHz to obtain an independent measurement of the con- 
tinuum emission. 

We observed the nearby source J1849+670 (distance 
to target ^^3.9°) every 15 min for amplitude and phase 
calibration. The strong calibrators J1927+739, 3C273 
and 3C454.3 were used for bandpass calibration. Nep- 
tune and MWC349 were used for absolute flux calibra- 
tion. Pointing was checked every hour using the optical 
and radio methods (jCorder et al.|[2010l ). The data were 
edited and calibrated using a combination of MIRIAD 
and CASA tasks. The first and last edge channels where 
the bandpass deteriorated were flagged accordingly. The 
large primary beam of CARMA at the observing fre- 
quency (1') allowed us to partially overlap the two point- 
ings resulting in a 20% increase in sensitivity. We sub- 
tracted the continuum emission from our data in the 
uu— plane by fitting the continuum with a polynomial 
of order 1, using line- free channels to estimate the con- 
tinuum level. This is particularly important for the 
H1821-I-643 line emission, for which strong continuum 
from the jet is expected (see below). 

The naturally weighted visibilities were imaged and de- 
convolved using CLEAN to a residual of approximately 
la, where a is the noise level, in a box centered in 
H1821-I-643. The final images reach a sensitivity of 
Icr w 0.7 mJy per 120 km s^^ channel in each field, with 
a clean beam size of 2.5" x 2.1" and a position angle of 
95°. 

3. RESULTS 

3.1. CO and continuum emission in H1821+643 

Figure [T] shows the CO velocity-integrated emission 
map and the corresponding spectrum averaged over the 



emitting region towards H1821+643. A Gaussian fit to 
the spectrum indicates an emission line with a full-width 
at half maximum (FWHM) of 580 ± 140 km s'^ and 
centered a t 120 zfc 60 km s~^ fr om the optical redshift 
z = 0.297 (jSchneider et al.lll992D . By fitting an elliptical 
Gaussian profile to the integrated CO map, we find that 
the emission appears to be extended along the major axis 
with a scale of 4.0" ± 1.2" or 18 ± 5 kpc, however it is un- 
resolved along the minor axis (Position Angle 132° ±8°). 
By fixing the extent of the emission to that value, we 
estimate an integrated CO emission of 2.3±0.5 Jy km 
8-^. The source is centered at a = 18''21'"57.43" and 
6 = 64°20'34".80, with a positional uncertainty from the 
fit of 0".8, in agreement with the expected positional 
uncertainty given by AO = 0.5 x SNR x 9^^^^ = 1.1", 
where SNR^ 5 corresponds to the signal to noise ratio 
and 6'bcam ^ 2.2" is the half-power beam width (HPBW). 

H1821-I-643 has an integrated continuum emission of 
10.9 ± 0.7 mJy and 10.3 ± 0.5 mJy at 88.9 GHz and 
93.1 GHz, respectively. The source is centered at a = 
18''21™57.207" and 6 = 64°20'36."22, with an uncer- 
tainty in the fit of 0.06". The emission is slightly resolved 
with a deconvolved size of~ 0.4" — 1.0" or 1.5 — 4.3 kpc. 
As shown in Fig. [TJ the continuum emission can be easily 
identified with the lo cation of the radio core of this QSO 
(|Blundell et al.l[T996l) . 

3.2. Previous reported CO observations of H1821+643 

Previous searches for CO in this object wi th the IRAM 
30m telescope failed dete c t any emission (lAlloin et al.l 
IT99llCo^es et al.ll20Tll ). iCombes eTall (|20lH l reports 
a 3(T limit of 1.7 Jy km s~^ in a 300 km channel. 
In order to directly compare this upper limit with our 
observations, we deconvolved the visibilities to a coarser 
velocity resolution of 300 km s~^, matching the one used 
by Combes et al., around the central velocity of the line 
profile. Integrating spatially over the emission area and 
over the velocity range covered, we find an integrated 
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Fig. 2.— Spectral energy distribution of H1821+643. The black 
crosses show the data points from radio to optical wavelengths. 
The green and blue dotted lines represent the dust and radio emis- 
sion models fitted to the data, with the solid green line showing the 
combined model. The red and blue solid lines show the redshifted 
radio loud and radio quiet QSO template SEDs from Elvis et al.l 
l |1994i ) and lPolletta et al.1 KOOTT ). respectively, scaled to match the 
optical emission. The magenta solid line shows the redshifted tem- 
plate SED of MrK231, a typical starburst to QSO transition object. 

flux density of ^ 0.72 Jy km s^^ over a 300 km s~^ 
channel, which is about 2.3 times deeper that Combes 
et al. observations. Hence, this explains why the CO 
emission line was undetected previously. 

3.3. Gas off the nuclei: a gas-rich/gas-poor galaxy 
merger or tail-like structure? 

The CO emission line from H1821-f643 is located 
2.0" (about 1 synthesized beam; ~ 8.8 kpc) south-east of 
the radio-core position, and 2.8" (~ 12.3 kpc) south-east 
of the optical center (host galaxy). Based on anal ysis o f 
PSF subtracted HST optical images, Fl ovd et al.l (|2004l ) 
find a nebulous optical structure. Floyd et al., however, 
point out that such optical feature may correspond to an 
artifact of the PSF subtraction procedure. 

The CO position coincides with the location of this 
optical structure, giving support to the reality of such 
feature. Indeed, the bulk of the molecular gas is located 
at only 0.52", well within the uncertainty in the posi- 
tion, from an optically faint peak to its north-east (Fig. 
[1]). Thus, this suggest that the CO emission is possi- 
bly associated with either an optically faint but gas-rich 
companion galaxy that is merging with the giant ellipti- 
cal host, or with an asymmetric tail- like structure within 
the host galaxy reminiscent of a previous merger event. 
Note, however, that the powerful AGN is hosted by a gi- 
ant elliptical galaxy, and large amounts of molecular gas 
(M(H2)> 5 X 10^ M0) such as that observed here (see be - 
low) are rare in early- type galaxies (jCrocker et al.|[20Tll ). 
Only 3 out of 56 early-type galaxies in th e Virgo clus- 
ter h ave molecular gas masses > 10^ M0 (jYoung et al.l 
1201 ID . being statistically unlikely that the molecular gas 
observed in H1821-f 643 is related to the host galaxy. Al- 
though it is still plausible that the gas could be related 
to a tail, relic of a past merger, the substantial amounts 
of molecular gas, its offset location with respect to the 



elliptical host and the existence of a starburst component 
in the optical-to-radio SED (see below) suggest that the 
CO emission actually comes from a different galaxy. If 
the molecular gas is indeed related to the host galaxy, it 
would imply a highly asymmetric gas distribution, possi- 
bly associated with a previous interaction or spiral arm, 
in an unusually gas-rich elliptical galaxy, likely related 
with starburst activity. Deep high-resolution IR observa- 
tions, transparent to obscuration, are necessary to decide 
between both possibilities. 

The CO emission is also well aligned with the radio- 
axis to the south-east suggesting a physical association 
between the radio-plume emission, the diffuse optical 
source and the large amounts of m olecular gas. Accord- 
ing to lBlundell fc Rawlingsl (|2001[ ). the radio jet appears 
to have been pointed within 7° within the last 0.3 Gyr. 
Thus, the probability that the CO emission is by chance 
located within 7° is ~ 2 x 7/360 or 3.9%, where the fac- 
tor 2 accounts for both radio lobes (north and south). 
Such probability is similar to the one found based on the 
alignment between the radio-axis and the CO position 
angle. 

3.4. ISM properties 

The CO luminosity (in units of K km s ~^ pc^) can be 
derive d from the line intensity following ISolomon et al.l 
(|1997[) as, 

L'co - 3.25 X 10^ (1 + z)-3 Dl f Scodv (1) 

Jav 

where J^y Scodv is the velocity integrated CO flux den- 
sity in units of Jy km s~^, z/obs is the observing frequency 
in GHz, and Z^l is the luminosity distance in Mpc at the 
source redshift z. From our integrated CO flux, we derive 
L'^Q = (1.0 ± 0.2) X 10^° (K km s"! pc^) . The lumi- 
nosity of the CO 1 — emission line is commonly used 
to estimate the molecular gas mass, M(H2), by using the 
relation M(H2)=Xco-^co- ^co is the CO luminosity 
to gas mass conversion factor which we assume to be 
0.8 M0 (K km s~^ pc^)~^, as it was found f or local IR 
luminous galaxies (jDownes fc Solomonlll998D . We thus 
estimate a molecular gas mass of (8.0 ± 1.7) x 10^ M0 for 
the CO emitting source. This molecular gas mass is at 
the high end of the mass range found for local ULIRGs, 
and is typical of luminous submillir neter galaxies and 
massive disk galaxies at high-redshift (iGreve et al.ll2005l : 
ISolomon fc Vanden BoutI [20051 : IDaddi et al.ll2010D . 

3.5. Dynamical mass of the CO-bright source 

Due to the relatively low significance of the CO de- 
tection, it is not possible to determine the actual geom- 
etry of the CO emitting source. We estimate the dy- 
namical mass based on the CO profile using disk geom- 
etry, which is commonly assumed for CO sources. How- 
ever, we remark that our adoption of a disk geometry 
is not well established for our CO source and we use 
it only as a guide and for the purpose of comparison 
with other CO studies. In this case, the dynamical mass 
is given approximately by Mdyn ~ (R / G) {v^ot / sin(i))'^ . 
The source radius can be expressed in terms of its ma- 
jor axis, as i? ~ ainajor/2 = 9 kpc, and the rotational 
velocity is given by the CO line- width ^ 580 km s"""^. 
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Fig. 3. — Density of galaxies in the redshift slice z 0.25 — 0.35 in 
the H1821+643 field. The contour levels indicate the overdensity 
of galaxies with respect to the average number of galaxies found in 
a control field to the south. 

The CO emission source has a major axis at the redshift 
of the source of ~ 18 ± 5 kpc (at the Scr level) and a 
minor axis limited to the beam size 10 kpc). This 
allows us to put a limit on the inclination angle given by 
i = cos^^(aniinor/amajor) > 40°. This takcs into account 
the uncertainty in the major axis and the value for the 
inclination angle would only increase with smaller minor 
axis. With these values, we find a dynamical mass of 

< 1.7 X 10^2 Mq. 

The observed CO distribution is highly asymmetric 
with respect to the QSO position itself and actually 
barely includes it (see Figure 1). Based on the accu- 
rate o ptical photometry of H1821-I-643 from HST im- 
ages, (IFlovd et al.l 120041 ) estimates a total mass for the 
elliptical host galaxy of 2.1 x 10^^ Mq. Given the asym- 
metric distribution of CO-emitting gas close to the QSO 
host galaxy, it is indeed surprising that the dynamical 
mass range estimated for the CO source amounts up to 

< 80% of the mass of the elliptical host (-- 2.1 x 10^^ 
M©). As explained above, such large amounts of molec- 
ular gas are rare (but possible) in old giant ellipticals 
([Crocker et al.l 120111: lYoung et all 1201 ID . and thus the 
CO emission could potentially arise from an obscured, 
merging starburst galaxy, which in this case would have 
an important stellar component. On the other hand, it is 
interesting to note that if the CO profile would have been 
used to estimate the QSO host galaxy dynamical mass, 
we would have underestimated such value by at least 
20%. This is consistent with results from high-redshift 
disk galaxies, for which the CO-based dyn amical mass 
was f ound to be underestimated by ~ 30% (iDaddi et al.l 

IMol) . 

The potential discrepancy between Mdyn(CO) and 
Mqso-host in massive mergers offers a general caution- 
ary tale when it comes to CO-deduced dynamical mass 
esti mates of QSO hos t galaxies at the higher redshifts 
(e.g. IWang et aLllMol) . 

3.6. Spectral energy distribution 

Using radiative transfer models of starburst and AGN 
emission. iFarrah et al- (,2002.) studied the spectral energy 



distribution (SED) of H1821-I-643. These models give a 
good fit to the far-IR data points although they fail to de- 
scribe the emission from old stars at wavelengths shorter 
than 4 /im. Their best fit model to the IR emission of 
H1821-I-643 indicates a composite object, with a warm 
AGN component that provides about 60% of the total 
luminosity, and a starburst component with a cold dust 
temperature of 30 K and a dust mass of ~ 3 x 10^ 
Mq. Their model, however, does not make reference to 
the radio emission, which could contribute a significant 
amount to the emission at submillimeter wavelengths. 

The spatial coin cidence between the radio core 
(jBlundell et al.lll996il and the peak emission in our 93 
GHz continuum images of H1821-I-643 implies that most 
of this emission comes from the nuclei and that it is non- 
thermal. As shown in Figure [2J the fiux density at ~ 93 
GHz is compatible with the trend followed by the radio 
emission at lower frequencies up to 350 GHz, in support 
for the non-thermal nature of such emission. 

In order to account for the radio continuum emis- 
sion, a power-law parametrization was used, of the form 
Si, oc {v/vq)", while simultaneously modeling the far- 
IR emission with a multi-component gray-body spec- 
trum . Each dust component can be described by a 
model of the form S^ oc Bi,{Td)il - e"^-), where B^{Td) 
is the Planck function and Td is the dust tempera- 
ture. The dust optical depth, Tj^, is proportional to the 
dust mass, Md, and to the dust absorption coefficient 
k(i/) = [v/vqY ^ which we adopt to be equal to 0.04 cm^ 
gr~i at 250 GHz (iKruegel &: SiebenmorgenI 11994 ) with 
P = 2.0 (IPriddev k McMahonll200'lll 

The best fit model is shown in Figure [H A cold dust 
component is found with Td = 50 K and a dust mass of 
1.4 X 10* Mq. a second, warmer dust component was 
obtained with a dust temperature of 130 K, and a dust 
mass of ~ 1.0 x 10^ M©. Both components alone do not 
allow us to fit the shorter wavelengths data points, A = 
6 — 10 /im , and thus we added a third, hot component 
with a dust temperature of 430 K, which we expect to be 
biased by the unknown amount of emission from stars at 
these wavelengths. 

The radio emission is well described by a power-law 
with a slope of a —0.05 and a flux density of 16 mJy 
at 1.0 GHz. The dust m odel obtained is consistent with 
that found bv lFarrah et al. (2002), however the inclusion 
of the radio emission implied a change in the resulting 
cold gas temperature from 30 K to 50 K. 

The cold dust component emission is identified with 
vigorous starburst activity, possibly related to the CO 
emitting source, and the warm/hot dust component 
emission with AGN activity. We find a total IR luminos- 
ity Ls-ioooAim = 1-0 x 10^^ Lq, with the AGN contribut- 
ing with 5.3 x 10^^ Lq and the starburst with 4.7 x 10^^ 
Lq. Note that although the flat spectrum radio emis- 
sion contributes with about 1% of the total IR luminos- 
ity, it amounts to ^ 70% of the emission at 350 GHz. 
On the other hand, the starburst emission is negligible 
at 90 GHz, ^ 25 /iJy, consistent with our observations. 
Finally, note that the dust mass for the cold dust compo- 
nent and our CO measurements imply a gas to dust mass 
ratio of ~ 60, comparable to the valu es found for a sam- 
ple of local far-IR s elected galaxies (IDunne et al.ll2000l : 
ISeaauist et al.|[2004f ) and distant submillimeter galaxies 
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Fig. 4. — Mosaic of the 93 GHz continuum emission obtained by 
combining the data obtained in the two overlapping target fields. 
The flux density scale in terms of the central noise level (0.3 mjy) 
is shown to the right. The cyan circle represents the position of 
H1821+643. Yellow squares show the position of the galaxies with 
redshifts matching the frequency coverage of the CO line in the 
cluster {z = 0.2925 — 0.3015). Yellow and white circles represent 
galaxies with photometric redshifts in the range z = 0.25 — 0.35 (24 
galaxies) and outside that range (25 galaxies), respectively. The 
size of circles and boxes is 2 X 6beam ■ 

(iGreve et al.l2005tlKovacs et al.ll2006l:lMichatowski et al.l 
l2010f ), bu t somewhat low co mpared to values found for 
ULIRGs (jWilson et aLllMish . 

4. DISCUSSION 

4.1. Another possible example of jet-induced galaxy 
formation 

Simulations predict that galaxies receive a large frac- 
tion of their gas through cold flows of intergalactic gas. 
Such neutral gas is found in filaments and cools ef- 
fectively in massive galaxy clusters onto giant ellipti- 
cal galaxies at their centers (Fabian 1994). Feedback 
from the central radio source may provide a mechanism 
for reheating the gas. It has been proposed that the 
cooled gas in the inner regions of the cluster may con- 
verge toward the central black- hole or collapse to form 
molecular gas clou ds and stars (jPizzolato fc Sokeill2005l : 
iRussell et al.ll20iol ). CO observations of the central mas- 
sive galaxies in clusters indicate that the deposition of 
gas toward the c entral regions is indeed produced in - 
short timescales (lEdgel 120011 : iSalome fc Combes! 120031: 
lEdee fc Fravedl2003fl . 

It is well-known that II1821-I- 643 lies at the center 
of a massive cluster of galaxies (jLacv et al.l 1992) and 
therefore the cold accretion scenario in this case is possi- 
ble. The gas funneled to the center would feed the mas- 
sive black- hole, and an accretion disk insta bility would 
ignite the jet emission (jRussell et al.l l20l"ol) . However, 
to start the outburst a galaxy merger scenario is more 
plausible since it would explain the precessing jet axis, 
which is likely prod uced by a binary black-hole system 
(jRussell et al.ll20l"o[ ). Our observations suggest that the 
actual answer is more likely a combination of both hy- 
potheses. The CO emission is offset from the QSO nuclei 
and coincident with the position of a faint optical feature 
to the south-east and well-aligned (position and position 



angle) with the southern jet axis. 

X-ray observations indicate that the cluster around 
H1821-h643 has a cooling rate of 300 ± 100 Mq yr-\ 
whereas the QSO acc retion rate, 40 Mfr^ yr~^, is at half 
the Eddington hmit (jRussell et al.ll2010t ). This implies 
that ~ 3 X 10^^ Mq has cooled down in the typical 
timescale of 1 Gyr of this cluster and is potentially avail- 
able for molecular gas formation and material accretion 
onto the central black-hole. At this accretion rate, the 
QSO would only keep growing for a fairly short period of 
time, < 0.1 Gyr, since in this time it would grow larger 
than the upper limit that has been observed in the local 
Universe (jRussell et al.ll20l"of ). This puts a limit on the 
period in which the QSO is active and might have formed 
jets, and is consistent with the timescale of 0.3 Gyr in 
which the QSO jet appea rs to have remained within a 
7° cone (BlundcU & Ra wlingsll2(30ll) . This implies that 
the jet axis has been pointed in fairly the same direction 
during most of its existence. 

The significant alignment between the jet axis and 
the companion galaxy as well as the jet lifetime (the 
time in which the jet has been impacting the merg- 
ing galaxy) suggest that the jet could have a role 
in boo sting star formation within the CO emitting 
source (jKlamer et al.ll2004l : lFeain et al.ll2007l : lElbaz et all 
|2009( ). Strong jet-ISM interactions inducing highly- 
excited CO lines over large molecular gas reservoirs has 
been recent ly noted in 3C293. a pow erful FR II ra- 
dio galaxy (jPapadopoulos et aLl l2008bl) . though its ef- 
fect o n the SFR (very low in 3C 293) remains far from 
clear (jPapadopoulos etaP l2010( ). Examples of radio- 
loud QSOs where the jet aligns with the location of the 
bulk of molecular gas have been foun d at high-redshift 
(jKlamer et al.ll200i : Hlbaz et al.ll2009t) . 

4.2. Searching for cold flows of gas as the origin for 
CO abundance 

In order to test how the gas is cooling onto the massive 
galaxies in this galaxy cluster, we mapped the distribu- 
tion of galaxies, which should roughly trace the filamen- 
tary structure and the path of intergalactic gas into the 
cluster core. We selected galaxies from the SDSS with 
photometric redshift redshifts between 0.25 and 0.35 and 
brighter than i' = 23.0 in an area of ^ 5.2' x 5.2' centered 
in II182 1-1-643. A density map from the selected galax- 
ies was constructed by creating a grid of cells and com- 
puting the number density as the 10th nearest neighbor 
galaxy to each cell position. Figure |3| shows the density 
map highlighting the position of the CO detected ob- 
jects. This map is given in terms of the average density, 
Sq = 0.81 arcmin"^, found in a control field to the south 
of the H1821+643 field, and thus represents the degree 
of overdensity {6/6o) with respect to this average. 

H1821-I-643 is well within the densest part of the clus- 
ter. Figure 131 suggests that there is density sub-structure 
associated with H1821-I-643, with the stronger density 
peaks following across the south-west, south-east and 
north-east of this source. This can also be seen in the 
distribution of galaxies shown in Fig. H) 

To investigate whether substantial amounts of gas have 
been deposited and condensed in molecular form onto 
other galaxies in the cluster, we performed a stacking 
analysis of the CO emission using the galaxies that have 
accurate spectroscopic redshift measurements. For this. 
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Fig. 5. — Left: Combined CO emission integrated over 400 km s^^ from the six galaxies that have a spectroscopic redshift measurement. 
Contours represent the fiux density in steps of la = 0.27 mjy beam~^ starting at ±2(7. Right: Combined 3.3 mm continuum emission 
over 23 galaxies with photometric redshifts in the range z = 0.25 — 0.35 within 30" from the center of the field. Contours represent the flux 
density in steps of Icr = 58 fijy bcam^-"^ starting at ±2(t. 



we selected galaxies that lie within a distance of 40" from 
of our target fields central pointing and in redshift range 
covered by the observing frequency of our observations, 
z = 0.2925 — 0.30 15, from the spectroscopic samples of 
I Tripp et al.l (|1998[ ) and references therein. With this cri- 
teria only six galaxies were selected, most of them lying in 
the outskirts of the image where the sensitivity decreases, 
as shown in Fig. We extracted 400 km s~^ single- 
channel maps centered at the position of each galaxy and 
with a central frequency corresponding to the galaxies 
redshift (e.g. i^co/(l + z)). We deconvolved the visibil- 
ities and the resulting dirty maps were thus stacked to 
compute a noise- weighted average map. We did not de- 
tect individual sources or their combined emission (Fig. 
4). We put a 3cr limit for their average emission of 0.81 
mJy per 400 km s~^, which constrains the gas mass con- 
tained on average in this galaxies to < 1.1 x 10^ M0. 
This is consistent with the molecular gas masses found 
in elliptical and lenticular galaxies in the local volume, 
which range between 2 x 10^ — 5 x 10^ Mq (jCrocker et al.l 
1201 It ), and excludes the possibility that gas- rich spirals, 
which typically have molecular gas masses of > 1 x 10^ 
Mq, are in the inner regions of the cluster. 

To check whether any possible bright continuum source 
could be contributing to the bright emission seen in FIR 
images, we take advantage of the comparatively high- 
resolution of the 3.3 mm images to search for continuum 
emission from the galaxies in this cluster field by us- 
ing stacking. To obtain better statistics, we increased 
our sample of galaxies in the cluster by selecting sources 
from the SDSS photometric catalog with redshifts in the 
range z = 0.25 — 0.35, optical emission i^g < 23 and 
located within 30" from the central position of each tar- 
get field. This resulted in a sample of 23 galaxies. We 
combined all line free bands making a total bandwidth of 
4.125 GHz at ^ 91 GHz. Similar to the CO maps, we ex- 
tracted continuum maps and deconvolved the visibilities 
and combined the dirty maps, which were thus stacked 
in a noise weighted fashion. Again no significant detec- 
tion was found in the average image down to a 3cr limit 



of 175 /iJy. This implies that the only strong continuum 
emitter in the field is indeed H1821+643. 

5. CONCLUSIONS 

We have presented the detection of ^^CO J = 1 — 
line emission towards the template ULIRG-to-QSO tran- 
sition object [H89]182H-643. The extended CO emis- 
sion is not associated with the QSO or host galaxies 
but is offset by about 8.8 kpc to the south-east. The 
molecular gas reservoir is estimated to be 8 x 10^ 
Mq. The CO emission coincides with the position of an 
optical feature suggested by HST optical imaging. We 
argue that H1821-I-643 could either represent a merger 
between a gas-poor, giant elliptical and a gas-rich com- 
panion galaxy or that the CO emission could arise from 
a tail-like gas-rich structure within the giant elliptical 
galaxy. 

Based on the optical emission, H1821-j-643 was found 
to have a total mass of ^ 2.1 x 10"'^^ M©. However, we 
find that the dynamical mass found for the CO source 
(offset from the elliptical host) based on its line profile 
can amount up to 80% of the mass of the elliptical host. 
On one hand, this suggest that the CO source should be 
a different galaxy that the elliptical one because other- 
wise its large mass would imply an highly asymmetric 
stellar distribution for the elliptical host. On the other 
hand, this results suggest that if we had measured the 
elliptical dynamical mass based on the CO profile, as- 
suming it is well centered in the elliptical, we would have 
underestimated its value. 

We revised the dual nature of the SED of H182H-643 
(starburst and QSO), and thus associate a cold dust 
(~ 50 K), massive starburst component with the CO 
emitting source, and the warm dust (~130 K), AGN com- 
ponent with the nuclei and host galaxy. We find that the 
inclusion of a flat radio spectrum component, important 
even up to 350 GHz, implies a higher temperature for 
the cold dust component respect to previous estimates. 

Based on the projected density map of galaxies in the 
surrounding cluster field, we find that the field shows an 
overdensity peak to the south-east of the H182H-643 sys- 
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tem, arguing in favor of cold flows of gas toward the inner 
regions of the cluster along this path. This hypothesis 
needs to be confirmed with sensitive CO and HI obser- 
vations, in order to reveal the nature of the CO emitting 
source and to directly measure the total amount of neu- 
tral gas deposited in the central regions of the cluster 
and QSO. 

The stacking technique was used to estimate the com- 
bined CO emission from galaxies in the field with spec- 
troscopic redshifts. This allowed us to put an upper limit 
to the molecular gas mass of 1.1 x 10^ Mq, comparable 
to that of typical elliptical galaxies. Additionally, the 
stacking analysis was used to search for 93 GHz contin- 
uum emission from a sample of photometrically-selected 
sources with redshifts in the range between 0.25 — 0.35. 
We found no strong continuum emitters down to a 3a 
limit of 175 fi Jy. The lack of significant amounts of gas 
in cluster field galaxies hints at the importance of in- 
teractions in the central object that can trigger the gas 
condensation and star formation. 

The results exposed above, in which the bulk of the 
molecular gas is well ofi^ the central QSO exemplify the 



importance of performing a revision of the local ULIRG- 
to-QSO transition objects with high spatial resolution 
observations. Many of these objects may have been mis- 
classified as post merger systems; and perhaps more im- 
portantly they are ideal targets to investigate unrelaxed 
gas-star configurations in which dynamical masses may 
have been underestimated. 
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